A relation between the dark mass of elliptical galaxies and their shape 
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It has been reported that there appears to be httle non-baryonic dark matter in some elliptical 
galaxies. We have studied a large number of elliptical galaxies and found a clear, previously unknown, 
correlation between their dark matter content and the ellipticity of their visible shape. At equal 
luminosities, flattened medium-size elliptical galaxies have on the average five times larger total 
mass than round ones. Such a correlation provides a new testing ground for models of dark matter 
and galaxy formation. 



Dark matter cosmology has been successful at deter- 
mining many large-scale features of the universe \V\ but 
open questions remain at the galactic and semi-galactic 
scales [2, Q. There, galaxies are described as made of 
large dark matter halos surrounding smaller luminous 
components. A search for relationships between the ob- 
served luminous matter of a galaxy and its dark mat- 
ter content can potentially advance our understanding of 
dark matter cosmology at this scale. Useful and intrigu- 
ing empirical correlations have been found, for example 
the TuUy-Fisher relation [j] for spiral galaxies, which 
relates their rotation speed to their absolute luminos- 
ity, or for elliptical galaxies, the Fundamental Plane |i5l] 
which includes the Faber- Jackson [d| and Kormendy [7] 
relations and links their effective radius Reff, their sur- 
face brightness and a, the statistical dispersion of stellar 
velocities. However, none of these involve directly the 
dark matter content of galaxies. There is also the puz- 
zling observation of little dark matter content in some 
elliptical galaxies [8|. We report here on a search for a 
correlation between the most obvious visible character- 
istic of elliptical galaxies, the ellipticity, on which their 
Hubble classification is solely based, and their relative 
amount of dark matter, expressed as the galactic total 
mass (dark-|- luminous) normalized to luminosity [M/L). 

Although elliptical galaxies can have tri-axial ellip- 
soidal shapes and ellipticities depending on radii, to first 
approximation their geometric shapes can generally be 
simply parameterized as oblate ellipsoids of constant el- 
lipticities (e). These span a large range, from spheres 
(e = 0) to flattened ellipsoids (e — 0.7), providing us 
with a group of smoothly varying shapes. However, only 
the projection of e on our observation plane (the appar- 
ent ellipticity, Sapp) can be measured, their true elliptic- 
ity {etrue) remaining elusive. In addition the M/L ratio 
for ellipticals is hard to measure accurately [8|. To over- 
come these difficulties, we select large and homogeneous 
samples of elliptical galaxies for which M/L has been 
extracted using different methods. Effects of galaxy pe- 
culiarities are minimized by a homogeneity requirement 
and suppressed statistically. The projection problem is 
addressed statistically. We select 255 different galax- 
ies for a total of 685 data points from 40 publications. 
Each provides M/L (or related ratios) for at least sev- 



eral galaxies, allowing us to study subgroups of data with 
consistent M/L extraction. From each publication, we 
select a sample of elliptical galaxies as homogeneous as 
possible to avoid diluting (by including effects particu- 
lar to a galaxy) or biasing any possible correlation (via 
systematic effects from a class of galaxies, e.g. compact 
or giant ellipticals). Only medium-size elliptical galax- 
ies are considered. These tend to have nearly isotropic 
random velocities (contrary to e.g. giant ellipticals). We 
require undisturbed galaxies to avoid galaxy-galaxy in- 
teraction from invalidating the method used to extract 
M/L (e.g. virial theorem, equations of hydrostatic equi- 
librium or strong lensing equations) . The strict selection 
applied to samples of local galaxies must be relaxed for 
samples of distant galaxies, since those are not as well 
characterized. (Distant galaxies are typically considered 
in strong lensing studies or in studies of the Fundamen- 
tal Plane time-evolution, see e.g. [9|). The local and 
distant rejection criteria are as follows (based on either 
NASA/IPAC Extragalactic Database (NED) 1^ or the 
publication from which M/L originates): 
Classes of rejected local galaxies: 

•Lenticular galaxies (SO- type). 

•Transition- type (E-I-), spiral galaxies. 

•Galaxies having an Active Galactic Nucleus (AGN) 
since, for mature galaxies, it may reveal that the galaxies 
have been disturbed recently. 

•Seyfert galaxies and BL Lacertae objects, for the same 
reasons as AGN. 

•Peculiar galaxies and galaxies listed in the Arp Atlas 
of Peculiar Galaxies. 

•Galaxies with HH emission, since the presence of HH 
regions is peculiar for ellipticals. Furthermore, it may 
reveal a recent disturbance. It may also bias AI/L since 
newly formed blue stars increase the luminosity. 

•Galaxies with low-ionization nuclear emission-line re- 
gions (LINER) since they may be due to AGN or star 
births, in both cases a sign that the elliptical galaxies 
have been disturbed. 

•Compact ellipticals since they belong to a different 
class of ellipticals. 

•Supergiant ellipticals (cD), giant ellipticals (D), 
Brightest Cluster Galaxies (BrClG) since they belong to 
different classes of ellipticals. Furthermore, giant ellipti- 



cals are characterized by more anisotropic random veloci- 
ties and tend to be triaxial. In addition, the dark matter 
of the cluster or group to which these galaxies belong 
may contribute to the galaxy dark matter content [11 ^ . 

•E? galaxies. We assume that this lack of nomencla- 
ture knowledge reflects poor measurements and may con- 
taminate the sample with non-ellipticals. 

•Bright X-ray and very faint X-ray galaxies. Very faint 
X-ray galaxies show signs of disturbed hydrostatic equi- 
librium ll|. In the case of bright X-ray galaxies, the hot 
gas property merges with the extragalactic gas of the 
group/cluster, indicating some influence from the envi- 
ronment and a contribution of the group/cluster dark 
matter to the galaxy M/L ll| . In addition, these galax- 
ies tend to be cD, D or BrClG. 

We keep low excitation radio galaxies and weak 
emission-line radio-galaxies as we have seen no reason to 
reject them. In some cases, when the statistics are low, 
we have relaxed some of these selection criteria (typically 
keeping LINER galaxies). In that case, the analysis is 
considered to be of lower reliability. This, together with 
the smaller statistics, makes these analyses contribute 
little to the final result. 
Classes of rejected distant galaxies: 

•Massive galaxies, typically with Af > 5 x IO^-'^Mq, to 
reduce the amount of cD, D or BrClG galaxies. 

•Galaxies with velocity dispersions a < 225 km.s^^, if 
no distinction is made in the publications between SO and 
elliptical galaxies, or if the classification is not reliable 
enough. This is to suppress possible SO contamination, 
since SO tend to have a < 225 krn.s"^. 

We remark that the rejection criteria have been chosen 
before performing the analysis. In this sense, the analysis 
is blind and as free of subjective bias as possible. 

We use data from six different categories of methods 
used to deduce M/L. This mitigates the possibility of a 
systematic methodological bias. These are: 

•Virial theorem: We employ virial data from 8 publi- 
cations |9J,ll^|l3|,ll5J. For cases where published analyses 
do not allow for ellipticity, we correct them using the ten- 
sor expression given in |13| . (This important correction 
derived analytically was experimentally verified in [14j.) 

•Stellar dynamics modeling: We use 11 publications 
la . llTj . An advantage of this method is that Stme can 
be inferred and used directly in calculations. 

•Interstellar gas X-ray emission data: The 2 publica- 



tions we used |11| assume spherical symmetry. We can 
partly account for ellipticity by replacing Newton's Shell 
Theorem (used in the derivation of the equation of hy- 
drostatic equilibrium) by a relation including ellipticity. 
This yields a correction numerically similar to the tensor 
correction [13j to the scalar virial theorem. 

•Planetary Nebulae and Globular Clusters data: The 
4 publications used here [8|, |l8| assume spherical symme- 
try. To account for ellipticity, we note that the ellipticity 
corrections to the virial and Shell theorems are close and 



apply the same correction. 

•Gas disk dynamics: We use 3 publications [19ll20l|. 

•Strong lensing: We use 16 publications [2, l2l|, 123 ]. 
Effects of £ are expected to be small 2l| and no correction 
is applied. 

Clearly, from such varied methods of determining 
M/L, different correlations with ellipticity can result 
depending on the treatment of e, systematic modeling 
biases, or sampling biases. For example, strong lens- 
ing tends to sample more massive galaxies, whose cor- 
responding M/L{e) correlation strength may differ from 
those of medium-mass galaxies. It is to minimize such 
particular biases that we employ all available methods. 

Linear fits of M/L to e are carried out for each of 40 
homogeneous samples. The fit results are then combined, 
which reveals a significant correlation between M/L and 
e. Results from 4 out of the 40 samples are shown on 
FIG. [H The large scatter results from the random pro- 
jection of the elliptical galaxies on our observation plane. 
However, all slopes indicate a positive correlation, with 
4.4, 4.3, 1.3 and 1.2 ct deviation from zero, from top left 
to bottom right panels respectively. (Fits of M/L with- 
out accounting for e, i.e. using straight horizontal lines, 
yield x^ /ndf values that are 1.5, 1.3, 1.1 and 1.1 time 
larger, respectively.) 
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FIG. 1. Apparent £ 

Galactic mass over luminosity in solar M0/L0 units versus 
ellipticity from four publications using different methods. We 
show some of the most statistically precise and systematically 
accurate results within each method. Top left: virial theorem 
[la |. Top right: strong lensing Q] (we plot 4 times the total 
mass over the stellar mass AMtot/M^ — M/L). Bottom Left: 
stellar orbit modeling \VT\. Bottom right: gas disk dynamics 



A difficulty with elliptical galaxies is that we observe 
only their projection, with their true ellipticities being 
unknown. However, the orientation of galaxies with re- 
spect to us is random. Hence, the Stme distribution can 
be computed from the Sapp distribution and the projec- 
tion corrected on a statistical basis. This correction is 
independent of the M/L extraction method if the ellip- 
ticity distributions in the galaxy samples are the same. 
We assume so and apply an identical correction to all 
the data sets (except for those already using etme, as for 
example in the bottom right plot of FIG. [1]). We use the 
£app distribution and M/L versus Sapp correlation ob- 
tained from [l3| to estimate the correction because it is 
the second largest sample of galaxies (the largest sample 
is from [15| but the authors indicate it may be slightly bi- 
ased). The relation between Strue and Sapp for an oblate 
spheroid viewed at an angle i is: 



i-va 



etrue)^sin^i -I- cos^i. 



(1) 



We assume a Gaussian distribution for etme- Its charac- 
teristics are chosen so that the simulated Sapp distribution 
reproduces the observed one. The projection correction 
depends on the relation between M/L and etme- The 
simplest assumption is a linear relation. We can also 
use a Bose-Einstein function to avoid possible unphysi- 
cal negative M/L near etme — 0. The distribution of 
M/L with Sapp can then be obtained from the Gaus- 
sian distribution transformed using Eq. ([1]) in which i 
is random. We also add a random shift toward smaller 
£app values to simulate the rounding effect of detector 
resolution. The initial linear or Bose-Einstein function 
values have been chosen so that the distribution of M/L 
with Sapp reproduces the experimental results. The fit- 
ted slope compared to the initial function's slope provides 
the projection correction. It increases ^^ \ the slope 
of M/L with e, by a factor 1.9±0.3. We must apply this 
factor to M/L ratios extracted using Sapp when they are 
combined to ratios extracted using Strue- 

Combining all results requires some caution. Most 
but not all M/L are extracted in the blue part of the 
electromagnetic spectrum, the B-band. There are scale 
factors between B-band results and those extracted in 
other frequency bands. Other scaling factors occur if 
authors used different Hubble constant values. Finally, 
sometimes Mtot/M^ is provided (M* is the stellar mass) 
rather than M/L (see e.g. the top right plot of FIG. [J) 
but Mtot/Mif is approximately proportional to M/L in 
most studies. To address these points, we normalize to 
M/L{eapp = 0.3) = 8M0/L0 = AMtot/M^{eapp = 0.3). 
The values M/L = SMq/Lq and M^/L = 4M0/LQ are 
typical in the B-band for elliptical galaxies. We normalize 
at eapp = 0.3 because the Sapp distribution peaks there 
and it avoids uncertain extrapolations which would be 
needed if, for example, we were to normalize at Eapp = 0. 
Another caution point is that some publications use sim- 
ilar methods to extract M/L and these have quoted re- 



sults for common galaxies. For example, the galaxy 
NGC7619 is used in 4 of the 8 publications employing 
the virial theorem. Consequently, their results may be 
correlated and bias the average. In addition some of the 
M/L extractions are, in our context, more reliable than 
others. To address this, the uncertainties are multiplied 
by a reliability factor. (This prescription allows some 
subjectivity but happens to have small influence, the less 
reliable results having usually also less statistics.) For 
extractions using similar methods, we count the com- 
mon galaxies and increase each uncertainty accordingly 
assuming that they are statistically dominated. (In do- 
ing so, we must account for the reliability factors since 
results are weighted by them when combined.) This pro- 
cedure assumes that results using a given method and an 
identical galaxy sample are perfectly correlated, which is 
not exact since within similar methods different assump- 
tions, proflles, input distributions, etc... are used. The 
speciflc analyses also differ in details and the data qual- 
ity varies as the publications used span a range of 30 
years. Thus, our procedure overestimates the uncertain- 
ties. This is partly mended by fitting 



d(M/L) 
de 



versus ra- 
dius with a one parameter function (i.e. ^ is taken 
to be constant with radius) and scaling the uncertainties 
so that x^ /ndf — 1. Forcing x^ /ndf = 1 is justified be- 
cause the results are from either different methods (for 
which we can assume that their systematic uncertainties 
are uncorrelated) or, for results using the same method, 
because we have accounted for the correlation. Although 
accounting for common galaxies and different reliability 
is necessary to average as accurately as possible, the net 
effect happens to be numerically unimportant. Finally, 
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can be 
,. This 



for M/L extracted using strong lensing, 
computed using either Sapp or mass ellipticity e 
later is the ellipticity of the modeled total mass distribu- 
tion. We use Eapp, as we do for the other five methods. 
We note however that 8 out of the 10 strong lensing pub- 
lications for which Smass is available display a stronger 
M/L versus Smass correlation, i.e. L < 



dCa 



dCn 



The slopes y are shown in FIG. [5J plotted for 
clarity against the approximate average radii at which 
the M/L values have been extracted. (Some papers pro- 
vide M/L calculated at different radii, so several points 
may correspond to one publication. Also, a few outlying 
points are outside the graph range. Consequently, there 
are more than 40 points in FIG. [2J) One notices sys- 
tematic shifts between results using different methods, 
most noticeably between the strong lensing and virial re- 
sults. This may be due to the different radii at which 
the M/L values are extracted. Also, a looser selection 
is applied on the distant -and thus less characterized- 
strong lensing galaxies, so contamination may dilute the 
correlation. In addition, no ellipticity correction is ap- 
plied because their extracted M/L are less sensitive to e. 
Interestingly, if we were to consider — — '—^ rather than 
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Slopes -2 — ^ versus approximate radii (units of Reff) at 
which the M/L are extracted. The six symbols distinguish 
the methods used to obtain M/L. The dash hue indicates 
the weighted average value after accounting for correlations. 



— — — , the strong lensing and the corrected virial results 
would be consistent. Finally, strong lensing galaxies tend 
to be more massive, so the comparison with our samples 
of local galaxies is not direct. Nonetheless, all the aver- 
aged \ from the six different methods are non-zero 
and positive. They average to, 

(^T^) = (17-42 ±4.30)Mo/Lo, 

a statistically meaningful positive signal, indicating that 
M/L and e correlate. We note that, without the selec- 
tion criteria, the otherwise clear correlation is obscured 
by large fluctuations, which may explain why it has not 
been noticed earlier. We verified that it disappears when 
no selection is applied. The correlation scales with our 
normalization MjLieapp — 0.3) = SMq/L©. Compared 
to SMq/Lq, the slope is large. 

An MjL correlation with e might reveal a significant 
bias in the data and/or methods. However we have made 
a thorough investigation of the various inter-dependences 
of the many variables characterizing elliptical galaxies 
and found that they could not explain the correlation 
seen. In addition, the six methods that have been used 
to extract MjL are independent. This indicates that 
the correlation is physical rather than a methodological, 
observational or instrumental bias. Its large magnitude 
makes it a feature that future descriptions of elliptical 
galaxies -including their dark matter content- must be 
able to reproduce. Numerically, with our normalization 
M/L{eapp = 0.3) = 8M0/L0, the correlation indicates 



that the most flattened medium-size elliptical galaxies 
have M/L averaging to (15.O±1.8)M0/L0 and the round 
ones to (2.8±1.2)M0/L0. Since in the B-band, the stel- 
lar mass over luminosity is M^/L ~ 4M0/L0, this im- 
plies that within the radius investigated (typically 0.1 to 
1 Reff), round medium-size elliptical galaxies have small 
amounts of dark matter (including in particular those 
considered in [8| ) . This is puzzling in light of the conven- 
tional model of cosmological structure formation, which 
requires dark matter seeds to trigger galaxy formations. 
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